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Abstract 

Abstract 

In this paper, the possibility to construct a path integral formalism by 
using the Hubbard operators as field dynamical variables is investigated. By 
means of arguments coming from the Faddeev-Jackiw symplectic Lagrangian 
formalism as well as from the Hamiltonian Dirac method, it can be shown that 
it is not possible to define a classical dynamics consistent with the full algebra 
of the Hubbard X-operators. Moreover, from the Faddeev-Jackiw symplec- 
tic algorithm, and in order to satisfy the Hubbard X-operators commutation 
rules, it is possible to determine the number of constraint that must be in- 
cluded in a classical dynamical model. Following this approach it remains 
clear how the constraint conditions that must be introduced in the classical 
Lagrangian formulation, are weaker than the constraint conditions imposed 
by the full Hubbard operators algebra. The consequence of this fact is ana- 
lyzed in the context of the path integral formalism. Finally, in the framework 
of the perturbative theory, the diagrammatic and the Feynman rules of the 



model are discussed. 
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I. INTRODUCTION 



The Hubbard X-operators M are suitable to give a powerfull framework in which the 
elementary excitations in solids can be explained. The use of X-operators is also relevant 
when electronic correlations are taken into account. This is the scenery in which High-T c 
superconductivity occurs, and so the main reason why the Hubbard operators algebra is so 
interesting at the present time. 

The algebra of the Hubbard X-operators is completely defined by: 

a) the commutation rules 

[Xf , X] s ] = SijiS^Xf - 8 a5 xf) (1.1) 

b) the completness condition 

xf + + xr = i (i.2) 

c) the multiplication rules for a given site 

XfTf = 5^Xf . (1.3) 

From now on and for simplicity we consider the case in which the indices a,/3 can only 
take the values + and — , and so the Hubbard X-operators are boson-like operators of the 
SU(2) algebra. The spin s = 1/2 is naturally contained in this case. 

It is easy to show that the equations (1.3) are not all independent, and so the full 
information contained in the algebra can be recovered from the equations (1.1), (1-2) and 
the following three independent equations 

X~+X ++ - X~ + = (1.4a) 
X+'X" -X+" = (1.4b) 

X + ~X- + -X ++ = (1.4c) 

Consequently, the full algebra given by equations (1.1)-(1.3) is equivalent to the commu- 
tation rules (1.1), the completness condition (1.2) and the three conditions (1.4). 
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A many body theory constructed by using the Hubbard operators as field variables, 
requires the application of techniques used in quantum field theories. From this point of 
view it is necessary to formulate the Wick theorem for the case in which the field operators 
are neither usual fermion nor bosons. Progress in this directions where done ||, but the 
problem is still open. 

Like in quantum field theories, another way to attach the problem is via the path inte- 
gral formulation. It is important to say that a suitable path integral formulation must be 
independent of a given representation. On the other hand it must be written in terms of an 
effective action with a well defined dynamics. This last point of view will be adopted in the 
present paper. 

The paper is organized as follows. In section II and III, by using the Faddeev-Jackiw (FJ) 
Lagrangian method 0, a general treatment for first-order Lagrangian systems containing 
the Hubbard operators as dynamical variables is given. A family of Lagrangian describing 
these dynamical systems is found. The use of these classical Lagrangians in a path integral 
quantization formalism is also analyzed. Strong arguments can be given showing that it is 
not possible to include the full Hubbard algebra (1.1)-(1.3) in a classical dynamical model. 
In section IV, we confront our results with others previously given in the literature. In 
section V, the diagrammatic and the Feynman rules for the model are constructed. Finally, 
conclussions and discussions are given in section VI. 

II. CLASSICAL LAGRANGIAN AND DYNAMICAL MODEL 

One of the traditional approaches to study the quantization of spin systems or t-J model 
in which the Hubbard operator algebra takes place, is to consider the constrained systems 
from the point of view of coherent state phase path integration. Also the usual Dirac's 
Hamiltonian method for constrained systems by considering slave boson or fermion repre- 
sentation is frequently used. 

By writing a family of first-order classical Lagrangian directly in terms of the four Hub- 
bard operators, our main purpose is to obtain information about the kind and the number of 
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constraints present in these models. In this way it is possible to obtain a response about how 
many information contained in the algebra (1.1)-(1.3) can be introduced at the classical level. 
This approach requires the introduction of a suitable set of constraints, a priori unknown, 
that must be determined later on. To this purpose it is useful to use the FJ Lagrangian 
method Therefore, we briefly indroduce some definitions and key equations. 

As is well known, the FJ symplectic quantization method is formulated on actions only 
containing first-order time-derivatives. The most general first-order Lagrangian is specified 
in terms of two arbitrary functionals Ka{h a ) and V^(/x), and is given by 

L^ A ^ A )=f, A K A ^ A )-\^{^). (2.1) 

The functionals, K A (fi A ) are the components of the canonical one- form K(n) = 
K A (/j,)djj, A and the functional V^(/x) is the symplectic potential. The general compound 
index A runs in the different ranges of the complete set of variables that defines the extended 
configuration space. 

The Euler-Lagrange equations of motion obtained from (2.1) are: 

X>^ B -^- = o. (2.2) 

The elements of the symplectic matrix M A b(h) are the components of the sympletic 
two- form M(n) = dK(fi). The exterior derivative of the canonical one-form K(/i) is written 
as the generalized curl constructed with partial derivatives and so, the components are given 
by: 



dK B 8K a 



When the symplectic matrix M A b is non-singular, from the equations of motion (2.2) 
result 



As the symplectic potential is just the Hamiltonian of the system, the equation (2.4) is 
written 



/A V = (2.5) 



.^4 R 



9V(o) 



dfx 



B 



where 

(M AB ) _1 , (2.6) 



are the generalized brackets defined in the FJ symplectic formalism. 

It is easy to show that the elements (M AB ) _1 of the inverse of the symplectic matrix 
Mab correspond to the Dirac brackets [0 of the theory. 

Transition to the quantum theory is realized as usual replacing classical fields by quantum 
field operators acting on the Hilbert space, where quantum ordering and proper boundary 
conditions for the quantum field operators must be taken into account. Therefore, the 
predictions of both FJ and Dirac methods are equivalents. 

When the matrix M AB is singular, the constraints appear as algebraic relations and 
they are necessary to maintain the consistency of the field equations of motion. In such 
a case, there exist m (m < n) left (or right) zero-modes v a (a = l,...,m, A = l,...,n) of 
the supermatrix Mab, where each v a is a column vector with n + m entries v A . So the 
zero-modes verify the following equation 

Y,v A M AB = 0. (2.7) 

A 

From the equations of motion (2.2) we see that the quantities fl a are the true constraints 
in the FJ symplectic formalism, and they are given by 

n a = < A v «« = o . (2.8) 

Consequently, in a first iteration the constraints are written in the symplectic part of the 
Lagrangian by means of Lagrange multipliers as follows 

L« = ^a^)+e a ^a-V« , (2.9) 

where the new symplectic potential is by definition = V^°' ) | Q=0 . The partition jj A = 
{ip l , £ a ) and Ka = {a% , £l a ) has been made. So, the compound indices A, B run the set 
A = (i , a) and B — (j , b). 

6 



In each iterative procedure the configuration space is enlarged and the symplectic matrix 
is modified. When no new constraints are found the iterative procedure is finished. 

Now, we are going to apply the FJ quantization formalism to a dynamical model for the 
Hubbard operators. 

As it is well known in all the examples in which the field variables are the components 
of the spin operators, the starting point is to consider first-order Lagrangians. This also 
happens in the t-J model when it is written in slave boson or fermion representation |§. 
The FJ quantization algorithm is suitable to study this kind of dynamical systems described 
by constrained first-order Lagrangians in which the constraints play a crucial role. 

Therefore, in the case under consideration we assume that the first-order classical La- 
grangian as functional of the Hubbard operators is written as follows 



where H(X) is for instance the Hamiltonian of the Heisenberg model written in terms of the 
Hubbard operators. The site indices were dropped since they are irrelevant in the analysis 
we will develop. Without any difficulty the site indices can be opportunely included. 

In the equation (2.10) A a is an adequate set of Lagrange multipliers which allows the 
introduction of the constraints in the Lagrangian formalism. Q a (X) is a set of suitable 
unknown constraints, initially considered ad hoc in the Lagrangian. Both the constraints 
Q a (X) as well as the range of the index a must be determined by consistency. The coefficients 
a a p(X) = a*p a (X) are found in such a way that the algebra (1.1)-(1.3) for the Hubbard 
operators must be verified. 

Looking at equation (2.10) we see that the initial set of dynamical symplectic variables 
is defined by (X a/3 , A a ) and the symplectic potential is given by 



So, the symplectic matrix (2.3) obtained from the Lagrangian (2.10) is singular, therefore 
the constraints are obtained by using the equation (2.8) and they read 



L = a a p(X)X a P - H(X) - X a n a , 



(2.10) 



V (0) = H{X) + X a Q a . 



(2.11) 



0V(o) 



(2.12) 
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and the first-iterated Lagrangian writes 



LW=a a ^X)X^ + i a Q a -H(X). 



(2.13) 



The modified symplectic matrix associated to the Lagrangian (2.13) is 



M AB = 



( 



da up 



dX a P 



) 



(2.14) 







where the indices A = {(a/3), a} ,B = {(7<5),fe}. 

At this stage the problem is to determine which, and how many constraints can be 
deduced from the algorithm of the method in such a way to obtain a non-singular symplectic 
matrix. 

In this way, from the Lagrangian (2.13) the symplectic matric (2.14) is constructed and 
its inverse can be computed. By equating each elements (M AB ) _1 of the inverse of the 
symplectic matrix Mab to each one of the commutations rules (1.1), differential equations 
on the coefficients a a p{X) and on the constraints Vt a are obtained. 

As it can be seen, the dimension of the symplectic matrix (2.14) is 4 + a, where a 
enumerates the constraints. Because of the antisymmetric property of Mab the index a has 
even range. From the properties of this matrix we can conclude: 

i) If a > 4 or odd, the symplectic matrix is singular. 

ii) For a = 4 the symplectic matrix can be invertible, but it is not possible to obtain 
the commutation rules (1.1). The commutators obtained by using equation (2.6) vanish, 
independently of the value of the coefficients a a p{X). On the other hand, when the number 
of constraints equals the number of fields there is no dynamics. So, it is not possible by 
means of Lagrange multipliers to enforce the constraint (1.2) together with the other three 
conditions (1.4). 

Consequently, we must resign the introduction in a classical first-order Lagrangian of the 
complete information contained in the algebra (1.1)-(1.3). 

Then, the unique possibility is to have only two constraints. The equation (1.2) or com- 
pletness condition must be imposed accounting their physical meaning. It avoids at quantum 
level the configuration with doubly occupied sites. The remaining constraint can not be any 
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one of that given in (1.4), because the commutators (1.1) can not be recovered. Therefore, 
we can expect that the remaining constraint can be provided naturally by consistency, when 
the symplectic method is used. 

Consequently, we assume an arbitrary constraint Q = Q (X + ~ , X~ + , u), where u = 
X ++ — X . This assumption is not a restriction because, by the completness condition, 
the sum (X ++ + X ) is equal to one. From the requirement that the matrix elements 
of the inverse of the symplectic matrix (2.14) must be equal to each one of the Hubbard 
commutation rules (1.1), and by solving the differential equation on this constraint the 
solution we find is 

D, = X+-X-+ + i« 2 - 13 = , (2.15) 

where /3 is an arbitrary constant. 

We emphasize that the constraint (2.15) is not an imposition but appears naturally from 
our method. This is the unique possible constraint in order to satisfy the commutation 
rules and the completness condition. Of course in equation (2.15) there is less dynamical 
information than the contained in the three equations (1.4). 

We will discuss this point connected to the fact that the path integral for this kind of 
fields represents the system in some limit of the operatorial approach. 



III. DETERMINATION OF THE LAGRANGIAN COEFFICIENTS. 

The next step is to determine the functions a a/ 3(X) written in the Lagrangian (2.13). 
The two constraints Q a we must consider are given in equations (1.2) and (2.15). Once the 
symplectic matrix (2.14) is constructed its inverse can be computed. Taking into account 
the equation (2.6) and the commutation rules (1.1), by consistency the following differential 
equation is found, 

<9a+_ da u + _ r <9a_+ da u + <9a_+ <9a+ . 

2[ ^r~dx^ ]x - 2[ ^r"aFT ]x +[ dx^-dx^ ]u = t > (3 - 1} 

where a u = \{a ++ — a ). 
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We assume that the coefficients a + _ , a_ + and a u can be written as products of arbitrary 
functions of the u variable by polynomials in the X + ~ and X~ + variables . For simplicity 
we try to look for a particular family of solutions by taking first-order polynomials in the 
X + ~ and X~ + variables, i.e 

a+ _ = f(u) [e + bX + ~ + cX~ + ] , (3.2a) 
a_ + = a\_ = f*(u)[e* + c*X + ' + b*X~ + } , (3.2b) 

a u = h(u) [p + qX + - + rX~ + ] , (3.2c) 

where the constant coefficients p , q , r , e , b and c are arbitrary ones. 

Once the expressions (3.2) are introduced in the equation (3.1) by straightforward com- 
putation we find 

ph{u) = iph{u))* (3.3a) 
qh(u) = (rh(u))* (3.3b) 
qh(u) = e-j- (3.3c) 

tt Hi 

cf(u) - c*f*(u) = 2tlm cf = 2*^^ (3.3d) 

with the conditions 6 = 0, and being a an arbitrary integration constant. 

Consequently, the equations (3.2) for the Lagrangian coefficients and (3.3), determine a 
family of Lagrangians compatible with the commutation rules (1.1), the completness condi- 
tion (1.2) and the constraint (2.15). 

Not losing generality, in equations (3.2b) and (3.3d) we can choose c = i and the function 
f{u) results 

/(«) 



4/3 — m 2 
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and so two different families of solutions are obtained: 
i) If e = the solution reads 



a n = - a — ) = K u ) > (3.4c) 

where h(u) is an arbitrary real function which also can be taken equal to zero, 
ii) If e 7^ the solution reads 

a +- = Jr^ (1 + iX ~ +) • (3 - 5a) 

a u = i(a ++ - a__) = %)[1 + X + ~ + X" + ] , (3.5c) 

where in this second case h(u) verifies the equation (3.3c). 

Really, the two different families of solutions (3.4) and (3.5) take into account the ma- 
jority of the significant cases. 

Finally, we note that making the following linear transformation to real variables 

(Si , S 2 , S 3 ) 

X + - =Si + iS 2 , (3.6a) 

X~ + = Si - iS 2 , (3.6b) 

X ++ - X— = 2S 3 , (3.6c) 

and by defining the vectors a = ( as 1 , as 2 , a>s 3 ) , V = ( ds 1 , <9s 2 , <9s 3 ) and S = 
(Si , S 2 , S 3 ) where 

11 



as 1 = «h — + a — h , 



(3.7a) 



a S2 = i(a+_ - a_ + ) , 



(3.7b) 



a S3 



a ++ — a , 



(3.7c) 



the equation (3.1) can be written in a more simple way and it reads 



(V x a).S = 1 . 



(3.8) 



The form of the differential equation (3.8) is equal to that obtained in Refs.[9,10]. Then, 
the fact that the kinetic term can be written as a function of a vector field a which satisfies 
the equation (3.8) is recovered. It must be noted that the equation (3.8) is a good definition 
for a curl on a S 2 manifold. Then, the equation (3.8) together with (2.15) written in terms 
of the new variables Si , S 2 and S3, allows us to write the kinetic term in the Lagrangian as 
the area of a sphere with /3 1 / 2 radious. This is the principal argument to say that (3 1 / 2 must 
be integer or half-integer. For a complete discussion about this argument see Refs.[9,10]. 

IV. A SIMPLE CASE AND ITS RELATION WITH PREVIOUS WORKS 

From section III, we can assert that a big family of Lagrangian really exists and any 
one of them can be considered as a good candidate for describing the dynamics contained 
in the commutation rules of the X-operators. The aim of this section is to discuss some 
important points by using explicitly one of the possible Lagrangians found in the previous 

section. Thus, by taking a ++ = a = in equation (3.4c) and calling a = —2s and j3 = s 2 , 

the Lagrangian (2.10 can be written 




(4.1) 



with the two constraints 




+■ 



+ 




(4.2a) 
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X ++ + X~~ = 1 . (4.2b) 

The equations (4.1) and (4.2) describe the classical dynamics of a system in which the 
commutation rules (1.1) are verified. 

We note that the same result also can be found by using the Dirac theory for constrained 
systems J7|. From this approach it is easy to show that the constraints given in equations 
(4.2) together with the constraints coming from the definition of the momentum of the X's 
variables, is a set of second class constraints. The Dirac brackets associated to this set of 
constraints are exactly the correct commutation rules for the Hubbard operators. 

Now, we are able to write the partition function by using the path integral Faddeev- 
Senjanovic approach [|TTJ and it reads 



Z = J DX 5[X + -X- + + i(X ++ - X—f - s 2 ] 5(X ++ + X— - 1) 

x expi [ dt L(X,X) , (4.3) 



where L(X,X) is given by (4.1). 

By integrating in the X variable we obtain for the partition function Z the following 
expression 

Z = J DX- + DX + -DX ++ 5[X + -X- + + ^(2X ++ - l) 2 ] 

: J dt L {1 \X,X) , (4.4) 



where 



x exp i 



L^X,X) = -i X+ - X ;;- X - +X ' + '-H { X). (4.5) 



Making in equation (4.5) the following change of variables, 



S, = X+ +X + , (4.6a) 



x+- - x-+ 

S 2 = — ^ , (4.6b) 



= \(2X+ + - 1) , (4.6c) 



13 



the functional integral (4.4) can be written as 



Z = J DS 5(S 2 + S 2 2 + Si - s 2 ) exp i J dt L^(S, S) , 



(4.7) 



where 



L^{S, S) 



S2S1 — S1S2 
s + S 3 



H(S) 



(4.8) 



where the constant Jacobian of the transformation (4.6) was absorbed in the functional 
integral measure. Therefore, the equation (4.7) for the partition function agrees with the 
expression (3.17) of Ref.[12], obtained by means of different arguments. 

Now, it is easy to show that this expression is consistent with the quantization of a spin 
system in the limit of large s. Applying again the Dirac theory, but now to the Lagrangian 
(4.8) with the constraints 



we find that the second class nature of the constraint defining Dirac brackets are again 
exactly the commutation rules (1.1) for the spin components. It is interesting to note that 
in the quantization procedure, the second class constraint (4.9) must be considered as a 
strong equation among operators. Then, 



From the comment given at the end of section III, it is known that the number s must 
be integer or half-integer. Consequently, in the equation (4.10) it is not possible to write s 2 



This is one of the important reasons which allows to ensure that in the path integral 
formalism for the spin systems, the information of large s approximation is contained from 
the begining. This fact is connected with our results making impossible the inclusion of 
the full X-operators algebra in a classical Lagrangian formalism, or equivalently in a path 
integral formulation. 



S | 2 = s 2 , 



(4.9) 



s 2 = s 2 i . 



(4.10) 



as s'(s + 1) with s integer or half-integer. 
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V. DIAGRAMMATIC AND FEYNMAN RULES 

Now, in order to obtain the diagrammatic and the Feynman rules for the model the 
perturbative treatment is analyzed. The starting point is to consider the following partition 
function 

Z = J DX + - DX-+ Du 5{X + -X- + + ^u 2 - {3) 

xexpi J dtL(X,X) , (5.1) 

where the integration on the (X ++ + X ) variable has been made by using the function 
5(X ++ + X— - 1). 

Consequently, the Lagrangian L(X,X) can be written: 

L = a + ^X + - + a„ + X~ + + a u ii - H(X + ~, X~ + , u) . (5.2) 

Taking into account the equations (3.4) for the coefficients, we consider the perturbative 

development for large value of the parameter f3. Therefore the non-polynomial Lagrangian 

(5.2), up to first order in reads: 

icy . . i 

L(X, X) = — (X~ + X+- - X + ~X-+) + a u ii + — u(X- + X+- - X + ~X~+) - H(X) . 

(5.3) 

In equation (5.3), we consider for the Hamiltonian H(X) the Heisenberg ferromagnetic 
form: 

H(X) = -\j(X + -X- + + X- + X + - + \uu) . (5.4) 

where J > 0. 

By using in the path integral (5.1) the Gaussian representation 

S(x) = lim — = exp( x 2 ) 

for the delta function, the partition function can be written in terms of an effective La- 
grangian as follows 
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Z = J DV expi £ dt L eff (V) . 



(5.5) 



In the equation (5.5), the effective Lagrangian L e ^ (V) is written in terms of an extended 
complex vector field V whose components are given by 



V = (X + ~ , X~ + , u) , 



and it can be partitioned as follows 



L*ff = L^(V) + L^(V) + L (4 \V) 



(5.6) 



As it is usual the quadratic part L^ 2 \V) of the effective Lagrangian defines the free prop- 
agator of the model, and the remaining parts L^(V) and L^{V) represent the interaction 
vertices, i.e the three and four legs vertices of the model respectively. So, from equation 
(5.5) it can be seen that the quantum problem remains defined in terms of a path integral 
which contains the three independent fields X + ~ , X~ + and u. 

In equation (5.6) the quadratic part L^{V) is given by, 



L {2 \V) = -V a D a ^, 



(5.7) 



where 



D, 



a/3 





■^9 t + l + J 



4/3 



d t + ! + J 








V 



(5.8) 





ad t + £ + iJ 

The simplest case in which a u = h(u) = au (where a is an arbitrary constant) has been 
considered when the matrix (5.8) was computed. 

The matrix D a p is Hermitian and non-degenerate, and so the propagador (Z)^)^ 1 in the 
[q, u] - space can be evaluated and it results, 

/ . Aa _. \ 



(D al3 ) 1 (uj,Uj') 



4/3 



V 














i 



(5.10) 
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We note that J q is the Fourier transforms of J^- = J only if i, j are nearest neighbor 



sites. 
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The three and four legs vertices are respectively given by the parts 

L^\V) = ^K Pl V a V^V\ (5.11) 



L^iV) = ±\ a/3lS V a V^V s , (5.12) 



where 



\ a fa = -uj)5(-uj + uj' +uj")5(-q + q +q") (5.13) 

K/3-yS = -^—(>{uj + uj + uj + uj )d{q + q + q + q ) a — (j — 'y — — 3, (5.14) 
2a 



\ 1 r / ' " "' \ r/ ' it in. . 

K/3-yS = o(— a; + a; + uj + uj )d{—q + q + q +q ) a = 1, (J = 2,^ = o = 3 

a 

and all the permutations, (5.15) 

\ 4 , / // ///. . / ii in. inn 

K/3-yS = 0{-uj + uj -uj +uj )0(—q + q - q +q J a = l,,a = 2,7 = l 

a 

and all the permutations. (5.16) 

From the above results we can see that for a = —2\f]5 = —2s and by choosing for the 
parameter a the value a — where z is the number of nearest neighbor sites, the matrix 
element reads 

(D 12 r = (T[X+-(c)X7 + (c')]) = - _ ^ - — 5{u - J)8{q - q) . (5.17) 

The above equation gives precisely the magnon propagator of the usual spin-wave theory. 

From (5.10), it can be seen that the longitudinal mode {T[u u}) has a pole on the imag- 
inary axis. This non-physical mode is related with the fact that there is no longitudinal 
dynamics in the lowest order of the spin-wave theory of the Heisenberg ferromagnetism. So, 
not losing physical information, also can be taken a u = 0. 

By computing the propagator and vertices for the solution (3.5) at the same perturbative 
order, it is easy to show that the same results are obtained. In particular the free propagator 
takes the form (5.8) with a = 0. 

In a further work under preparation, we will apply our perturbative approach to the 
renormalization and dumping of magnon energies. 
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VI. CONCLUSIONS AND DISCUSSIONS 



In this paper a new discussion about the path integral formalism for dynamical systems 
written in terms of Hubbard operators is done. 

If it could have been possible, this path integral must have contained the full algebra 
(1.1)-(1.3) for the X-operators. Using the Faddeev-Jackiw symplectic formalism we have 
shown that this proposal is not possible, and in order to satisfy the commutation relations 
rules (1.1) we have resigned to include the complete information contained in the X-operator 
algebra (1.1)-(1.3). 

By consistency of the formalism and in order to satisfy the Hubbard commutation rules 
we have found the number of constraint conditions. From our point of view and in a total 
independent way we arrive to a path integral which is consistent with those obtained by 
means of coherent states method. 

We have also shown that this path integral for the spin system case, is valid in the 
large spin s limit. Then, our conclussion assert that this limit is closely related with the 
impossibility to include the full algebra of the Hubbard X-operators in a classical dynamics. 

On the basis of our path integral formulation we present the diagrammatic and the 
Feynman rules for the perturbation theory. We have shown that our free theory is consistent 
with the results provided by the lowest order of the spin-wave theory. 

Finally, we can emphasize that from our approach a large family of kinetic terms of 
effective Lagrangians can be found. Some of them can be related with previous Lagrangians 
obtained by different methods. 
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